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ABSTRACT 

We present a study on the origin of the metallicity evolution of the intra-cluster 
medium (ICM) by applying a semi-analytic model of galaxy formation to iV-body/SPH 
(smoothed particle hydrodynamic) non-radiative numerical simulations of clusters of 
galaxies. The semi-analytic model includes gas cooling, star formation, supernovae 
feedback and metal enrichment, and is linked to the diffuse gas of the underlying simu- 
lations so that the chemical properties of gas particles are dynamically and consistently 
generated from stars in the galaxies. This hybrid model let us have information on the 
spatial distribution of metals in the ICM. The results obtained for a set of clusters 
with virial masses of ~ 1.5 x 10 15 h^ 1 M© contribute to the theoretical interpretation 
of recent observational X-ray data, which indicate a decrease of the average iron con- 
tent of the intra-cluster gas with increasing redshift. We find that this evolution arises 
mainly as a result of a progressive increase of the iron abundance within ~ 0.15i? v ir- 
The clusters have been considerably enriched by z ~ 1 with very low contribution 
from recent star formation. Low entropy gas that has been enriched at high redshift 
sinks to the cluster centre contributing to the evolution of the metallicity profiles. 

Key words: galaxies: clusters: general - methods: numerical - cosmology: observa- 
tions - X-rays: galaxies: clusters 



1 INTRODUCTION 

The intra-cluster medium (ICM) is a hot (~ 10 7 K) and dif- 
fuse gas contained within the deep potential well of clusters 
of galaxies, which constitute the largest virialized structures 
of the Universe. This gas radiates energy through thermal 
Bremsstrahlung, which is detected in the X-ray band of the 
electromagnetic spectrum, characterized by the presence of 
emission lines from highly ionized iron (Fe XXV and Fe 
XXVI) at 6.6 - 7 keV. These lines allow to obtain infor- 
mation about the metallicity of the ICM; the average Fe 
abundance is measured to be around 0.5Zq in the central 
regions of local X-ray clusters^ 

Supernovae explosions are the typical sources of met- 
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1 We note that this value is different from the one found in the 
literature, w hich is based on the photo sferic solar Fe abundance 
published in lAnders fe Grevessei ]l989h . We are using the more 
recent s olar value of iron abun dance by number (Fe/H)Q = 2.82 X 
10" 5 bv lAsplund et al.l ||2005|) . 



als that contaminate the ICM, being generated through 
subsequen t episodes of star formation in cluster galax- 
ies (e.g., iRenzinil 1 19971 ). although other sources may 
have a non negligible contribution to the metal bud- 
get of the ICM, like hypernovae assoc iated with popula- 
tion t ype HI stars (ILoewensteinl l200l|) an d intra-cluster 
stars (|Zaritskv, Gonzalez fe Zabludofj 20041 ). The metals 
produced enrich the ICM by different processes that af- 
fect the heavy eleme n t distribution, like galactic winds 
jHeckman et alj l200d. ISpringel fe Hernquistl 120031 ). ram 
pressu re stripping dKapferer et al.ll2007l ). and active galactic 
nuclei (|Moll et al.ll2006l ). The knowledge of the history of the 
chemical enrichment of the ICM contributes to our under- 
standing on the way and epoch in which the galaxy clusters 
have formed, the history of star formation in cluster galaxies 
and the associated supernovae rates, and the physical pro- 
cesses involved in the diffusion of energy and metals within 
the hot intra-cluster gas. Obse rvations of radial abundance 
profiles of different elements (|Pe Grandi fe Mol endi 2001; 
iTamura et"afll2004 IVikhlinin et alj|2005l ) provide valuable 
constraints on the physical processes involved in the chem- 
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ical enrichment of the ICM, especially those related to the 
feedback mechanisms that inject metals into the diffuse 
phase. 

A recent analysis of Chandra X-ray spectra of 56 clus- 
ters within the redshift range 0.3 < z < 1.3 spanning tem- 
peratures 3 < kT < 18 keV a llows to trace the ev olution of 
the iron content of the ICM l|Balestra et al.ll2007T ). The re- 
sults are based on the estimated average emission-weighted 
(EW) iron abundance of the ICM within the inner region 
of the cluster, delimited by ~ 15 — 30 per cent of the virial 
radius, R v i r . They find that, for z > 0.5, the mean abun- 
dance of the ICM is approximately constant, with a value of 
Zf c rs QAZq. However, at lower redshifts (0.3 < z < 0.5), 
the EW iron abundance is significantly higher, reaching a 
value of Zf ~ 0.64 Z@. The parametrization of these results 
implies that the iron abundance at the present epoch is a fac- 
tor of ~ 2 larger than at z ~ 1.2. The se results are supported 
by the more recent study made by iMaughan et al. (2008) 
based on a larger sample of galaxy clusters at 0.1 < z < 1.3 
observed with Chandra ACIS-I. 

From the current observed ICM metallicity evolution, 
it is not clear whether it arises as a result of a change in 
the iron content of the intracluster gas due to enhanced star 
formation, or as a consequence of redistribution of metals 
within the central regions; the later would be favoured by 
dynamical processes which transfer the chemically enriched 
gas from the intergalactic medium of the cluster galaxies to 
the hot intracluster gas. 

Since measurements of ICM abundances are emission- 
weighted, the mean abundance of the intracluster gas con- 
tained within a given radius is biased to higher values when 
considering the central regions of the clusters. This bias is 
expected to be higher in cool core clusters because of the 
peaks in their surface b rightness and iron abundan ce that 
characterize them (e.g., |Pe Grandi fc Molendil 12001! ). How- 
ever, the observed evolution of the ICM iron abundance 
cannot be entirely attributed to a possible decrease of the 
fract i on of cool-core cluster s at high redshift (|Balestra et al.l 
12001 IMaughan et alj|2003) . 

Some studies on the origin of the central peak in 
the metal distribution of clusters with a cool core sug- 
gest that it has been likely produced by th e brightest clus- 
ter galaxy after the cluster w as assembled l|Bohringer et al.l 
12004 |De Grandi et al.ll20o4 ). For this mechanism to be ef- 
fective, large enrichment times (4 — 10 Gyrs) are necessary. 
This conclusion has been achieved from the combination of 
observational results and population synthesis and chemical 
enrichment models. This possible mechanism for the forma- 
tion of the central peak has been complemented with the 
effect of the AGN-induced flows as a possible transporta- 
tion processes that have mixed the metal s into the ICM 
l|Rebusco et al l 120051 ; iRoediger etlll 120071 ). trying to ex- 
plain the observed difference between the broad ICM metal- 
licity profiles and the much narrower stellar light profile of 
the central galaxy. 

Different theoretical approaches have been used so far 
aimed at studying the processes of metal enrichment of the 
ICM. By using observed and modelled rates of type la su- 
pernoyae (SN e la) and core- collapsed supernovae (SNe CC), 
lEttoril (|2005T l has evaluated their contribution to the history 
of metal accumulation in the ICM, finding that they provide 
a good qualitative agreement with observations of the overall 



decline in abundance with redshift thanks to the large de- 
lay times for SNe la adopted. However, the total amount of 
iron is systematically lower than observed values in the red- 
shift range < z < 1. This model indicates that half of the 
iron mass observed locally is produced by SNe la, with SNe 
CC products becoming dominant at higher redshifts . Dif- 
ferent conclusions are obtained by iLoewensteinl (|2006l ) who 
uses a similar kind of model, implying instead that synthe- 
sis of cluster iron was dominated by SNe CC and/or SNIa 
with short delay times whose progenitors originated during 
a phase of rapid, top-heavy star formation. 

The impact on the ICM chemical evolution of envi- 
ronmental effects, such as ram pressure, and tidal or vis- 
cous stripping has also been analyzed. Calura, Matteucci 
& Tozzi (2007) use chemical evolution models applied to 
the different morphological types of cluster galaxies. They 
find that iron enriched gas ejected by galactic winds aris- 
ing from ellipticals accounts for the ICM iron abundances 
at z > 0.5, whereas the gas stripped from the disk galax- 
ies turning into SO through interaction of the interstellar 
mediun with the ICM accounts for the ICM metallicity at 
2 < 0.5. The joint contribution of these processes can ex- 
plain the increase i n the ICM abundance in the low redshift 
range observed by iBalestra et alj l|2007h . By using a com- 
plete different approach based on a combination of TV-body 
with a phenomenological galaxy formation model, which 
provides the ejected galactic matter that is then included 
into a hydrodynamic ICM simulation. iKapferer et all (120071 ) 
find that the mass loss triggered by ram pressure stripping is 
more dominant than galactic winds since a higher redshift 
[z < 2); however, their enrichment model cannot account 
for the observed evolution of the the mean ICM metallicity 
with redshift. 

From an hybrid model that combines an 7V-body/SPH 
non-radiative simulation of a gala xy clu s ter w ith a semi- 
analytic model of galaxy formation, ICoral l|2006h has shown 
that dynamical effects play a major role in the development 
of central iron abundance peaks. In the proposed scenario, 
the intergalactic medium of infalling cluster galaxies, that 
has been primarily enriched at high redshifts and then trans- 
fered to the ICM, is subsequently driven to the cluster centre 
by bulk motions in the intracluster gas. 

Using cosmolog ical hydrodynamical simulations, 
iTornatore et al.1 ((2007) also find that gas dynamical effects, 
related to gas mixing, galactic winds and ram pressure 
stripping, play an important role in determining the distri- 
bution and relative abundance of different chemical species, 
being the typical age of the ICM enrichment of z ~ 0.5. 
However, the influence of an excess of star formation at low 
redshift is more impor t ant t han the one obtained with the 
hybrid model of lCoral (|2006l ). which gives quite low recent 
star formation for cluster galaxies. 

For the theoretical interpretation of the observational 
results concerning the ICM metallicity evolution, we apply 
here the hybrid model of the chemical enrichment of the ICM 
implemented by Cora (2006). Although this hybrid approach 
does not couple metal productio n to the gas actually c ool- 
ing in radiative simulations fe.g.. ITornatore et al.|[2007l ). its 
novel feature is that the chemical properties of the diffuse gas 
in the underlying non-radiative TV-body/SPH cluster simu- 
lations are generated from metals ejected from the galaxies, 
consistent with the modelling of the semi-analytic model. 
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These metals are then carried around and mixed by the hy- 
dro-dynamic processes during cluster formation, allowing us 
to study the evolution of the spatial distribution of metals in 
the ICM. We compare the results of our hybrid model with 
th e mean iron abundanc es of the ICM recently determined 
bv lBalestra et ah 

I (|2007h . with the aim of investigating the 
connection between the development of the chemical abun- 
dance patterns that characterize the ICM and the observed 
average metallicity evolution. 

This paper is organised as follows. Section [5] briefly de- 
scribes the hybrid model used to study the ICM chemical 
enrichment, presenting the changes applied to the version 
described in lCoral l|2006h . mainly related to the distribution 
of metals among gas particles; it also summarises the prop- 
erties of the hydrodynamical simulations used. Section [3] 
contains the analysis of the evolution of the ICM chemical 
enrichment considering both radial abundance profiles (Sec- 
tion 13. If) and mean metallicities (Section I3.2|l of the main 
cluster progenitor at different redshifts. In Section[3]we sum- 
marize our conclusions. 



2 MODEL OF THE ICM CHEMICAL 
ENRICHMENT 

We consider an hybrid model for studying the chemical en- 
richment of the ICM, which consists of a combination of non- 
radiative cosmological 7V-body/SPH simulations that con- 
tain galaxy clusters and a semi-analytic model of galaxy 
formation. The important feature of this hybrid model is the 
link between semi-analytic model results and the chemical 
enrichment of the diffuse gas component of the underlying 
A^-body/SPH simulation. We pollute gas particles with met- 
als ejected from the galaxies that have been generated by 
the semi-analytic model. By enriching the gas particles lo- 
cally around galaxies, we can account for the spreading and 
mixing of metals by hydrodynamical processes, thereby ob- 
taining a model for the evolution of the spatial distribution 
of metals in the ICM. 



2.1 Semi-analytic model of galaxy formation 

The se mi-analytic model used here is based on previous 
works (jSpringel et al.ll200ll : |Pe Lucia et all 120041 ) . but was 
extended with a new chemical implementation that tracks 
the abundance of different species produced by different 
sources. We briefly describe here the modelling of the circu- 
lation of metals among the diffe rent baryoni c components. 
Details of this model are given in ICoral (|2006l ). The simplic- 
ity of the semi-analytic model has the advantage of reach- 
ing a larger dynamic range than fully self-consistent hydro- 
simulations, at a far smaller computational cost. In partic- 
ular, it allows us to explore more easily the range of pa- 
rameters that characterize appropriate chemical enrichment 
models. The long cooling time of the bulk of the gas in rich 
clusters justifies the assumption of a non-radiative gas in 
the ./V-body/SPH simulations. However, this time-scale be- 
comes smaller than the Hubble time in the core of the clus- 
ter, where the densities are higher. Here is where the semi- 
analytic model plays an important part, taking into account 
radiative cooling, star formation, and chemical enrichment 



Table 1. Number of high resolution dark matter particles, 
Nh r , of the non-radiative cosmological hydro-simulations, and the 
virial mass, M v - lr , temperature, T v - lT and radius, R vlr , of the most 
massive cluster in each simulation. The simulation name used in 
iDolag et al. I j2005l) is given between brackets. 



Cluster 


N hT 


M viT 










[r'Mg] 


[keV] 


[h- 1 Mpc] 


CI (g51) 


2219034 


1.3 x 10 lb 


7.8 


2.3 


C2 (gl) 


4937886 


1.5 x 10 15 


8.4 


2.1 


C3 (g8) 


5602561 


1.7 x 10 15 


8.1 


2.8 



and energetic feedback from galaxies by including the effect 
of supernovae explosions. 

Dark matter haloes and substructures that emerge 
in the simulation are tracked by the semi -analytic code 
and u sed to generate the galaxy population (jSpringel et al.l 
l200l| y In this subhalo scheme, the virial mass of a dark mat- 
ter substructure is defined simply as the sum of the mass of 
its particles, and the virial radius is estimated by assuming 
that the subhalo has an overdensity of 200 with respect to 
the critical density. Since each particle in the ./V-body/SPH 
simulations is split into dark matter and gas, the identifi- 
cation of dark matter haloes is based on the dark matter 
particles with their masses corrected for the cosmic baryon 
fraction. This is done to get haloes defined as close as possi- 
ble to pure dark matter runs, on which semi-analytic models 
are usually calibrated, but still allows to make use of the hy- 
drodynamic treatment for the distribution of metals ejected 
by galaxies. 

We assume that the hot gas always has a distribution 
that parallels that of the dark matter halo. Its mass is ini- 
tially given by the baryon fraction of the virial mass of the 
dark matter halo, and subsequently modified by gas cooling, 
star formation and feedback processes. Stars can contami- 
nate the cold and hot gas because of mass loss during their 
stellar evolution and metal ejection at the end of their lives. 
The hot gas has primordial abundances initially (76 per cent 
of hydrogen and 24 per cent of helium), but becomes chem- 
ically enriched as a result of the direct ejecta from stellar 
mass loss, and the transfer of contaminated cold gas to the 
hot phase that is reheated by supernovae explosions. This 
chemical enrichment has a strong influence on the amount 
of hot gas that can cool, since we are using m etal depen- 
dent cooling rates (ISutherland fc Dopital 0*9931 ) . This pro- 
cess in turn influences the star formation activity which is 
ultimately responsible for the chemical pollution. 

The chemical model implemented in our semi-analytic 
code considers mass losses from stars in different mass 
ranges. Massive stars give raise to core collapse supernovae, 
which include those of type Ib/c and II, being the later the 
most abun dant on e s. We adopt stellar yields as given from 
models of iMarigol (2001) for low- and intermediate-mass 
stars (0.8M Q < M < 5 - 8M ) and iPortinari et al.l (|l998f ) 
for quasi-massive (5M© < M < 8Mq) and massive stars 
(8Mq < M < 120 M ). Ejecta from type la supernovae 
are also included, considering the n ucleosynthesis p r escrip - 
tions from the updated model W7 bv llwamoto et ail f l999i. 
We us e the stellar lifetime given bv lPadovani fc Matteuccl 
(1993) to model the return time-scale of the ejecta from 
all sources considered, being specially relevant for the sin- 
gle stars in the low- and intermediate-mass range and for 
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the progenitors of SNe la, that are characterized by a 
long delay time from the formation of the progenitor to 
the supernova explosion. We adopt the single degenerate 
model to estimate t he SN e la rate, following the scheme of 
iGreggio fc Renzinil i|l983l ). where type la supernovae origi- 
nate in binary systems whose components have masses be- 
tween 0.8 an d 8 Mq . Calculations are based o n the formalism 
described in lLia. Portinari fc Carrar 

3 <|2002h . but assuming 
that core collapse supernovae, originate from single stars 
with masses larger than 8 Mq . 

The mass range of secondary stars in binary systems 
gives explosion times for SNe la comprised between ~ 2.9 x 
10 7 and ~ 1.4 x 10 10 yrs, with the SNe la rate reaching a 
maximun within ~ 0.1 — 0.7 Gyr for a single stellar popu- 
lation. Thus, the ejection rate of elements produced by SNe 
la reaches a maximum at z ~ 3 — 5, while the contribution 
of SNe CC to the ICM chemical enrichment peaks at higher 
redshifts (z ~ 5 — 7) since SNe CC rate closely follows the 
star formation rate as a result of the short lifetimes of the 
stars involved. The peaks in the mass ejection rates are fol- 
lowed by a strong decline at lower redshifts for both types of 
SNe, such tha t the ongoi ng chemical contamination is quite 
low at z = (ICoralbOOd 1 ). The relative contribution of SNe 
CC and SNe la to the iron content of the cluster is already 
one-third at z ~ 2, as given by the accumulated iron masses 
provided by both sources. 

Some aspects of the semi-analytic model have been 
modif ied with respect to the version described in ICoral 
(2006). We consider here a Salpeter IMF normalized be- 
tween 0.1 and 100 Mq. In order to obtain a mass-metallicity 
relation who se evolution with redshift approaches to the ob- 
served one (|Erb et al.1 120061 ). half of the mass ejected by 
galaxies because of mass loss through stellar winds or su- 
pernovae explosions are transfered directly to the hot phase 
instead of being first incorporated to the cold gas. In this 
way, the cold gas associated to each galaxy is more gradually 
contaminated, instead of already achieving at high redshifts 
the higher metallicitiy values that correspond to the present 
epoch. 

The free parameters of the model regulate the way in 
which gas cooling, star formation, supernovae feedback and 
galaxy mergers proceed, and determine the circulation of 
metals among the differ ent baryonic components. They have 
been carefully tuned in lCoral (120061 ) to satisfy numerous ob- 
servational constraints (Milky Way properties, the luminos- 
ity function, the Tully-Fisher, color-magnitude and mass- 
metallicity relations). The feedback efficiency parameter, 
that regulates the amount of cold gas that is reheated by 
core collapse supernovae, is now assigned a smaller value 
of 0.1 that allows to recover a better agreement of model 
results with the observed slope of the Tully-Fisher relation. 

2.2 Non-radiative AT-body/SPH cluster 
simulations 

The semi-analytic model is applied to cosmological hydro- 
dynamical simulations. We consider a set of three simulated 
galaxy clusters, having virial mass in the range ~ (1 — 
2) x 10 15 h~ x Mq. These clusters have been initially selected 
from a dark matter s imulation of a cosm ological box, having 
size of 479 h _1 Mpc (|Yoshida et aUl200lD . for a cosmological 
model fi m =0.3, ft A =0.7, H c = 70kms _1 Mpc" 1 , fi b = 0.039 



for the baryon density parameter, and as = 0.9 for the 
normalization of the power spectrum. The Lagrangian re- 
gions surrounding the three selected clusters have been re- 
simulated at higher mass res olution by applying the Zoomed 
Initial Condition technique dTormen et al.lll997f ). The mass 
resolution is the same for all simulations, with masses of 
dark matter and gas particle of md m = 1.13 x 10 9 ft -1 M© 
and m gas = 1.69 x 10 8 Mq, respectively. As for the force 
resolution, the Plummer-equivalent gravitational softening 
is fixed at e = 5 h' 1 Kpc in physical units at redshift z ^ 5, 
while it switches to comoving units at higher redshifts. Ta- 
ble 1 gives the number of high resolution dark matter parti- 
cles of each simulation, which is equal to the number of gas 
particles, and the virial properties of its most massive clus- 
ter. The simulations have been car ried out using the Tree- 
SPH GADGET-2 code (|Springe]||2005l) . GADGET- 2 is a parallel 
Tree+SPH code with fully adaptive time-stepping, which in- 
cludes an integrat ion scheme which explicitly conserves en- 
ergy and entropy l|Springel fc Hernquisj|2002l ). The simula- 
tions considered here include only non-radiative physics and 
the original formulation of artificial viscosity within SPH. 
General properties of these simulated clusters, such as ra- 
dial profiles of velocity dispersion of dark matter and gas 
particle s, and of gas tempe rature, density and entropy, are 
given in lDolag et ahl (|2005l ). 

2.3 Hybrid model of the ICM metal enrichment 

The key physical processes that lead to the evolution of the 
iron content of the ICM are the production of metals and 
their deposition in the gas. The former is controlled by the 
prescriptions included in the semi-analythic model, while the 
later depends on the way in which chemical elements gener- 
ated by galaxies in the semi-analytic model are distributed 
among gas particles in the ./V-body/SPH cluster simulation, 
which thus carry the information on the spatial distribution 
of chemical and thermodynamical properties of the ICM. 
The metal deposition results from a complex combination 
of different processes that includes powerful outflows from 
star forming galaxies, generated by supernovae explosions 
and winds from young massive stars, dynamical stripping 
and gas mixing. We model the deposition of metals in a 
simple way. For each snapshot of the simulation, we iden- 
tify a set of iVgas gas particles within a sphere characterized 
by an enrichment radius, -R G nrich, centred on each galaxy. 
Metals ejected by each galaxy are spread among the iVg as 
gas particles. Both the choice of -Ronrich and the procedure 
used to distribute the chemical element s have been modified 
with respect to the version presented in Cora (2006). In that 
work, metals were uniformly distributed within a radius of 
100 Kpc for all galaxies. 

In the present implementation, the mass of iron ejected 
by the galaxy is distributed among its corresponding set 
of iVg as gas particles according to the same kernel used for 
SPH calculations in the Af-body/SPH simulation, that has 
been obtained from a B-spline. Besides, we now consider a 
physically motivated model for determining 7? onr i c h based 
on a treatment of galactic winds (Bertone, Stoehr & White 
2005). It follows the expansion of supernova-driven super- 
bubbles around several hundred thousand galaxies formed 
in a region of space with diameter 52/i _1 Kpc and mean 
density close to the mean density of the universe. About 



Dynamical Origin of the ICM Metallicity Evolution 5 



half of the enclosed galaxies at z = are field galaxies, 
while the rest are in groups and poor clusters. They solve 
the equation of motion for a spherical astrophysical blast- 
wave to follow the evolution of winds after they escape the 
visible regions of galaxies. Although the history of each wind 
depends on the properties of the pare nt galaxy and on the 
environment where the wind expands, IBertone et ail (|2005T l 
find small scatter in the distribution of the shock radii for 
different parameter choices of the model at a given redshift 
and as a function of the stellar mass of the galaxy. This en- 
sures that mean quantities of global properties of the winds 
well represent the general trend of the whole populatio n. 

Considering the results of IBertone et al.l l|2005h . we 
adopt the dependence of the mean bubble radius for all 
galaxies blowing a wind at z = 3, as a function of the stellar 
mass, M*, of the host galaxies (see their Figure 3). This de- 
pendence varies according to the parameters chosen for the 
wind model, such that the bubble radius ranges from « 100 
to 400 h' 1 Kpc for galaxies with M* « lO n M . In order 
to explain the origin ot the observed evolution of the iron 
content of the ICM, we first have to reproduce the iron spa- 
tial distribution at z = 0. We find that the curve for which 
galaxies with stellar mass M+ ps 10 11 M© have a bubble ra- 
dius of 100 Kpc allows to recover radial iron abundance 
profiles in good agreement with observations, as shown in 
Section f3.1l We fit this curve with a parametrization given 
by 

-Rcnrich = -623.13 + 133.01 logAf* - 6.13 log 2 M*. (1) 

Using a larger bubble radius for massive galaxies, o r the 
same radius for all galaxies (fixed at 100 /i -1 Kpc as in ICoral 
2006) give place to radial iron abundance profiles flatter than 
observed. 

It is important to note that we are not considering the 
evolution of the bub ble radius given by the treatement of 
IBertone et al.l (|2005T l , but only the relationship between this 
radius with the stellar mass of the host galaxy at a given 
redshift. The subsequent spreading and mixing of metals 
is carried out by gas particles themselves through hydro- 
dynamical processes. The consistent implementation in our 
model of the equations governing the wind evolution and 
their impact on the chemical pollution of the ICM will be 
presented in a forthcoming work. 



3 EVOLUTION OF ICM IRON CONTENT 
3.1 Radial abundance profiles 

Observed iron a bundances at different clustercentric radi i 
for local clusters l|Tamura et al. | |2004 IVikhlinin et"al]|2005h 
constitute an important constraint to our model. We eval- 
uate its capability to represent the chemical enrichment of 
the ICM by comparing the radial iron abundance profiles 
of the ICM constructed from simulated results at redshift 
2 = with this kind of observational data. 

Radial abundance profiles are built by dividing the 
volume limited by the cluster virial radius into concentric 
spherical shells centred on the dominant galaxy of the clus- 
ter. These shells have equal width till 0.1i? v ir and are log- 
arithmically spaced beyond this radius. For each shell, we 
estimate the X-rs.y emission-weighted mean iron abundance 
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Figure 1. EW radial iron abundance profiles for the three sim- 
ulated clusters (thick solid, dott ed and dashed l i nes). Symbols 
represent observational data from I Vikhlinin et al. (2005). 

relative to hydrogen, using the emission measure of each gas 
particle defined as EM = J n p n c dV, where n p and n c are 
the proton and electron densities, respectively, and V is the 
volume associated to the gas particle. The iron abundance 
is then expressed in terms of the sol ar value, for which w e 
have adopted the recent calibration of lAsplund et al. (2005). 
These mean abundances are assigned to the mean radius of 
each shell. 

Figure [1] shows the EW radial iron abundance profile 
of the ICM at z = for the three simulated clusters (thick 
lines). They are compared with a sample of nearby, relaxed 
galax y clusters observed with Chandra (jVikhlinin et al.l 
2005), with a range of average temperatures similar to that 
of the simulated clusters. The comparison of the model 
abundances with observed values is shown out to R500, 
which is the radius enclosing a mean overdensity of 500 with 
respect to the critical density at the redshift of the cluster. 
The local EW profile shows a good agreement with the trend 
traced by observations, except in the innermost radial bins, 
being at odds with self-consistent hydrodynamical simula- 
tions that include gas cooling. The presence of dark matter 
substructures with associated high density and chemically 
enriched gas makes the EW abundance profiles quite noisy 
at distances larger than ~ 0.2 R500 because of the strong de- 
pendence of the X - ray emissivity on the gas density. These 
peaks are not present in mass-weighted (MW) abundance 
profiles, which show almost the same trend than EW ones 
but are much smoother. 

The evolution of the radial abundance profiles is fol- 
lowed by identifying the main progenitor of the cluster at 
different redshifts. Then, we take into account the gas par- 
ticles that reside within the virial radius of these structures. 
The main progenitor of a given halo at each redshift is se- 
lected among the set of progenitors of that halo and defined 
as the substructure with the largest number of dark mat- 
ter particles. Figure [2] shows the evolution with redshift of 
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Figure 2. Evolution with redshift of EW radial iron abundance 
profiles. They are obtained by averaging at each concentric shell 
the abundance distribution of the three clusters analysed. The 
error bars around the average radial profile at z = denote the 
1-a standard deviation of the mean. 

EW radial iron abundance profile of the ICM. Each curve 
corresponds to a different redshift and has been obtained 
by considering the gas particles within each radial bin for 
all the three clusters, thus giving an average abundance dis- 
tribution. The thick solid line corresponds to the present 
epoch, being the average of the individual radial profiles 
represented by the thick curves in Figure [1] 

From Figure [2] we can appreciate that the EW iron 
abundances become progressively lower with increasing red- 
shift, with the main change occurring at radii less than 
~ 0.2i?soo (~ 0.15 J?vir), where the profiles become flatter. 
At larger radii, the shape of the abundance distributions at 
different redshifts keep themselves quite close to the profile 
at z = 0. That is, the evolution of the radial abundance pro- 
files from z ~ 1 to 2 = is mainly produced by an increment 
of the iron content in the innermost regions of the clusters. 

3.2 Mean metallicity 

Recent results on the evolution of the iron content of the 
ICM based on the analysis of Cha ndra JT-ray spectra o f 
56 clusters of galaxies at 2 > 0.3 (|Balestra et al.ll2007l) . 
support a decrease of metallicity with redshift that can be 
parametrized by a power law of the form ~ (l+z)~ 1-26 . This 
observed evolution implies that the average iron content of 
the ICM at the present epoch is a factor of ~ 2 larger than 
at z ~ 1.2. Since this analysis is restricted to the central 
regions (R ~ 0.15 — 0.3 iivir) of clusters, the evolution of the 
iron content might be interpreted as a result of ocurrence 
of evolution in the gradients of the metal distribution. This 
hypothesis becomes possible taking into account our results 
on the evolution of the EW radial iron abundance profiles 
shown in Figure [2] 

In order to compare our results with those of 
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Figure 3. Evolution of the EW mean iron abundance of the 
main progenitors of the galaxy clusters analysed (thin lines). 
These mean values have been estimated using gas particles within 
0.15i? v i r . The mean values for all the three clusters are shown at 
some redshifts (filled diamonds; error bars denote la confidence 
level). The corresponding MW mean iron abundances estimated 
within 0.15i? v i r and i? v i r are also shown (empty diamonds and 
empty squares, respectively ). These result s are compared with 
the observational data from Ifjalestra et al ] j2007l ) (filled circles) 
and the corresponding parametrization of the metallicity evolu- 
tion with redshift (thick solid line). The observational data have 
been multiplied by a factor of 1.66 in order to r efer the abun- 
dances to the solar values of lAsplund et aL | [|2005h . 

iBalestra etakl (|2007l) . we estimate the mean EW iron abun- 
dance from the iron content of gas particles contained within 
0.15 R v ir from the centre of the cluster or of its progenitor, 
depending on the redshift considered. Figure shows this 
comparison for each of the clusters analyzed, together with 
the corresponding mean values for all the three clusters es- 
timated at some redshifts (filled diamonds). Observational 
data is represented by filled circles and the thick solid line 
denotes the corresponding parametrization of the metallic- 
ity evolution. Model results closely follow the trend denoted 
by this parametrization within the whole redshift range con- 
sidered. The innermost area contained within 0.15 J? v ir man- 
ifests the strongest changes in the slope of the radial abun- 
dance profiles (see Figure [2}. Thus, the decrease of the ob- 
served EW mean iron abundances with increasing redshift 
can be connected to the flattening of the central part of the 
abundance profiles. 

The rate of change of the EW mean metallicity with 
redshift has a mild dependence with the region of the clus- 
ter considered, manifesting clear evolution even for the case 
where all gas within R v i r is taken into account (not shown 
in the figure). The low density and poorly contaminated 
gas particles lying at distances greater than ~ 0.15 iivir do 
not contribute much to the EW mean metallicity because 
of their low Jf-ray emissivity. Hence, EW mean abundances 
mainly provide information of the central parts of the clus- 
ters. 
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Figure 4. Evolution with redshift of iron mass ejected by galax- 
ies within 0.15 -R v j r (thin lines) and -R v i r (thick lines) below a 
given redshift for cluster CI, expressed in terms of the iron mass 
contained within those radii at z = 0. Different types of lines 
indicate the joint or separate contribution of supernovae type la 
and CC. 



Instead, estimations of MW mean iron abundances give 
equal weight to both high and low enriched gas particles 
since their masses do not change in a non-radiative SPH 
simulation. Therefore, MW mean iron abundances strongly 
depend on the region of the cluster considered. Taking into 
account gas particles within either 0.15 Tivir or 7? v i r leads to 
MW mean metallicity evolution with smaller values when 
all gas contained within R V [ T is involved in the calculations. 
This is a consequence of the presence of progressively less 
contaminated gas particles with increasing clustercentric ra- 
dius. Results corresponding to R v i r and 0.15 7? v ir are shown 
in Figure [3] by empty squares and empty diamonds, respec- 
tively. 

As a result of the evolution of the slope of the radial 
abundance profiles within 0.15i? v ir, the MW mean metal- 
licity do show stronger evolution with redshift only when 
these inner parts of the cluster are considered, becoming 
quite close to the behaviour given by EW values, as shown 
in Figure [3] Hence, the observed metallicity evolution can 
be explained as being produced by the combination of two 
effects. On one hand, the evolution on MW mean iron abun- 
dances estimated from gas contained within 7? v j r indicates 
that there is an overall increase of the iron content within 
the virial radius of the main progenitor of the cluster. This 
is also evident from the lower values displayed by the whole 
radial abundance profiles as redshift increases (see Figure [2]). 
On top of this, we find that the evolution of the slope of the 
central part of the iron abundance profiles (< 0.15 7? v ir) is 
the key ingredient that contributes to achieve the observed 
evolutionary trend for z < 1.4. 



3.3 Understanding the evolution of the ICM 
metallicity 

The joint evaluation of observational and simulated results 
supports the fact that the observed evolution of the ICM iron 
content till redshift z ~ 1.2 is mainly due to the evolution 
of the metallicity gradients of the intracluster gas character- 
ized by an increase of central abundances with decreasing 
redshift. Both the star formation rate and the mass accre- 
tion history of the cluste r may give r aise to such evolution. 
The study carried out bv lCoral l|2006l ). based on the analysis 
of a single simulated cluster, favours the scenario in which 
dynamical effects play the major role in the increment of the 
central abundances of the ICM, since the amount of star for- 
mation at z < 0.5 is not enough to account for the observed 
evolution of the iron content. 

Figure U shows the history of metal ejection since z < 
1.4. For the simulated cluster CI, we estimate the accumu- 
lated iron mass ejected into the ICM below a given redshift 
by galaxies contained within 0.05 i? v ir and R v i T ; it is ex- 
pressed in terms of the total mass of iron contained in the 
ICM within those radii at z = 0. We see that, for z < 1.4, 
the accumulated iron mass ejected from galaxies within R V [ T 
is ~ 10 per cent, while this fraction is reduced to ~ 4 — 5 
per cent when only galaxies within 0.15_R v ir are taken into 
account. In this last case, the ejected iron mass drops quite 
abruptly at z ~ 0.3, making negligible contribution to the 
iron budget of the ICM within the innermost regions. We 
also show separately the contribution of different types of 
supernovae. We find that ~ 90 per cent of the iron mass 
ejected since z < 1.4 is provided by SNe la. 

These values of mass of iron ejected can be considered as 
lower limits because of the low star formation rate of clus- 
ter galaxies at low redshifts in our model. This behaviour 
arises as a result of the suppression of gas cooling in the 
central galaxies of halos with K-ir > 350km s^ 1 and in satel- 
lite galaxies; they simply form stars till their reservoir of 
cold gas is exhausted. The first condition has been imposed 
to avoid too massive cluster central galaxies that are not 
consistent with observations. The reduction or supression 
of the co oling flow could be naturally produced by AGN 
feedback (ICroton et all 120061 . iBower et all 120061 ) which is 
not currently included in our model. As a result of the dif- 
ferent modelling, it is not surprising to find that the iron 
ejected within R v i T since z < 0.5 is a factor of two smaller 
than the es timations obt a ined by the completely different 
approach of ICalura et alj t|2007l ). based on chemical evolu- 
tion models and morphological information, or those aris- 
ing fro m self-consistent cosmo logical hydrodynamical simu- 
lations (jTornatore et al.ll2007l ). which has an excess of low- 
redshift star formation. 

The fact that, at low redshift, the ICM in our model 
receives a very small contribution of iron ejected by cluster 
galaxies, allows us to isolate the role of dynamical processes 
in the evolution of the ICM abundance patterns. With the 
aim of understanding the way in which gas dynamics con- 
tributes to develop the iron abundance profiles, we select 
gas particles located within a sphere of radius 0.15 i? v ir cen- 
tred on the cluster at z = 0. They are tracked back to 
high redshifts in order to analyse the evolution of the re- 
lationship between their iron content and entropy values. 
Figure [5] shows the dependence of such quantities for clus- 
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Figure 5. Evolution with redshift of the dependence between the entropy and iron abundance of gas particles that have distances from 
the cluster centre within a sphere of radius 0.15 -R v j r at z = 0. Results correspond to cluster CI. Redshifts are indicated in each panel. Gas 
particles are colour-coded according to their distances with respect to the centre of the main progenitor of the cluster at each redshift: 
within 0.02 ij v i r (red), in the ranges 0.02 — 0.05i? v i r (magenta), 0.05 — 0.15i? v i r (blue), 0.15 — 0.3i? v i r (cyan) and 0.3 — 1 -R v j r (green). 
Gas particles outside il v i r are showed in Figure [6] 
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Figure 6. Evolution with redshift of the dependence between the entropy and iron abundance of gas particles that are outside the virial 
radius of the main progenitor of the cluster and end up within a sphere of radius 0.15i? v i r centred at the cluster CI at z = 0. Redshifts 
are indicated in each panel. 



ter CI at redshifts z = 0, 0.5 and 1 for those gas particles 
that remain within the main progenitor of the cluster as the 
redshift increases. Dots in the entropy-metallicity planes are 
colour-coded according to the distance of gas particles with 
respect to the centre of the main progenitor of the cluster: 
wihin 0.02i? v ir (red), in the ranges 0.02 — 0.05 _R v ir (ma- 
genta), 0.05 - 0.15i?vir (blue), 0.15 - 0.3i? v ir (cyan) and 
0.3 — 1 i?vir (green). The behaviuor of gas particles lying 
outside the virial radius of the main progenitor is presented 
in Figure [6] at redshifts z = 0.1, 0.5 and 1. 

The distribution in the entropy-metallicity plane of gas 
particles located at clustercentric distances < 0.15_R v ir at 
2 = (right panel of Figure [5| simply reflecs the fact that 
the closer are the particles from the cluster centre, the lower 
are their en tropy values, acco rding to the entropy profile of 
the cluster (|Dolag et al.ll2005l ). Only ~ 19 per cent of these 
particles resides at distances < 0.15 iivir at z « 1 (left panel 
of Figure [5}, while a larger percentage (~ 58 per cent) is 
found outside the main progenitor of the cluster (left panel of 



Figure |SJ). This behaviour is consistent with the flattening of 
metallicity profiles as redshift increases (see Figure [2}. Both 
sets of gas particles, within and outside the virial radius 
of the main progenitor, have already achieved a high level 
of enrichment by z ~ 1. Around 14 — 15 per cent of gas 
associated to other progenitor halos is characterized by low 
entropy values (< lOkeVcm 2 ) which are not so common 
among gas particles lying within the main progenitor. 

As the cluster assembles, gas particles outside the main 
progenitor of the cluster are incorporated into it following 
the infalling substructures, which are stripped in the pro- 
cess. At z ~ 0.5, the percentage of gas particles contained 
within the virial radius of the main progenitor increases to 
~ 63 per cent, with ~ 40 lying at < 0.15 _R v ir- Very few par- 
ticles are found outside the main progenitor at z » 0.1, as 
can be seen from the right panel of Figure [6] A small frac- 
tion of gas is stripped quite early from the halo progenitors; 
they are then incorporated onto the cluster by some diffusive 
accretion mechanism, increasing their entropy as a result of 
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the shocks suffered during the process (|DoIag et al.ll2005l ). 
As time passes by, the enriched low entropy gas converges 
to the inner regions of the cluster. The entropy profile of 
the ICM regulates this contribution. Although we have non- 
radiative simulations, the temperature profile of the ICM at 
2 = de creases for radius sma ller than 0.05 Rvir. As dis- 
cussed bv lChurazov et al l (|2003l ). the presence of this dense 
and cool core only allows the penetration of the low entropy 
gas attached to the infalling dark matter substructures that 
has been enriched at high redshift when the iron ejecta from 
galaxies were larger. The mass accreted onto the main pro- 
genitor at a given time may vary among the different clus- 
ters analysed because of their particular formation history; 
however, the general features of the evolution in the entropy- 
metallicity plane is the same for all of them. 

From this analysis, we find that much of the central, 
most highly enriched gas originates in infalling substructures 
that are not the main cluster progenitor. This scenario of the 
origin of the iron abundance pattern is supported by the re- 
sults of a test case in which gas particles are contaminated 
only by galaxies contained within the main progenitor of 
the cluster. In such a restricted enrichment scheme, the iron 
abundances are very low and the shape of the profile be- 
comes quite flat in comparison with the one obtained when 
the contribution of all galaxies is taken into account. The 
dynamical process that involves the sinking of enriched low 
entropy gas explains the evolution suffered by the metallic- 
ity profiles in their central parts, which cannot be accounted 
for by recent enhanced star formation activity. We can see 
the steepening of the profiles as a consequence of the spatial 
redistribution of the intracluster gas with different level of 
chemical enrichment. 



4 CONCLUSIONS 

We have applied an hybrid model for metal enrichment of 
the ICM that combines hydrodynamical cosmological sim- 
ulations of galaxy clusters and a semi-analytic model of 
galaxy formation. This hybrid model has the special fea- 
ture of linking the metal production of galaxies in the semi- 
analytic model with the chemical enrichment of gas particles 
in the Af-body/SPH simulations. This allows us to follow the 
evolution with time of the spatial distribution of metals in 
the intracluster gas. 

We have adopted a scheme for metal spreading among 
gas particles which is based on the dynamics of the ex- 
pansion of supernova-driven s uperbubbles generated around 
galaxies l|Bertone et alj|2005l ). We find that this physically 
motivated estimation of the enrichment radius satisfies ob- 
served ICM metallicity properties both locally and at high 
redshifts, being in good agreement with iron abundance pro- 
files at z = and the evolution of the mean abundances 
with lookback time. The simulated mean metallicities within 
0.15 -Rvi r are in good agreeme nt with the observations anal- 
ysed by iBalestra et all |2007) in the whole redshift range 
considered, closely following the parametrization of observed 
data till z < 1.4. The results of this model let us conclude 
that: 

(i) The observed mean metallicity evolution 
l|Balestra et al.l 120071 ) can be explained as being pro- 
duced by the combination of an overall increase of the iron 



content within the virial radius of the main progenitor of 
the cluster, and the evolution of the slope of the central part 
of the iron abundance profiles (< 0.15 -Rvir). This result has 
been obtained from the analysis of both MW and EM iron 
abundance profiles and mean metallicities, with EW mean 
values mainly providing information of the central parts of 
the clusters. 

(ii) The metallicities of gas particles that at z = end 
up within a region delimited by a radius of ~ 0.15 -Rvir have 
been already achieved between z ~ 1 and ~ 0.5, when al- 
most half of the gas has not yet been accreted onto the main 
progenitor of the cluster, and the metal ejecta from galaxies 
was considerably higher with respect to the present epoch. 

(iii) As time passes by, the enriched low entropy gas, 
mainly attached to infalling substructures, sinks into the 
cluster and is mixed there contributing to develop the iron 
abundance profiles. The presence of a dense and cool core al- 
lows the penetration of this low entropy gas l|Churazov et al.l 
2003). Hence, the very central metallicity (< 0.02 -Rvir) orig- 
inates almost exclusively from gas in infalling substructures 
that converge to the cluster core. This dynamical process 
explains the evolution suffered by the metallicity profiles in 
their central parts, which cannot be accounted for by recent 
enhanced star formation activity. 

The shape of entropy profiles determine the develop- 
ment of the abundance profiles. The forecoming conclusions 
have been obtained for non-radiative simulations of galaxy 
clusters. The entropy profile ch ange substantia l ly for radia- 
tive runs, as has been shown bv lBorgani et al. (2005), with 
the main change occurring within ~ 0.1 -Rvir. For simula- 
tions that include gas cooling, star formation and feedback 
regulated by the velocity of galactic winds, the entropy pro- 
file monotonically decreases to the cluster centre, while in 
non-radiative simulations an almost isentropic core develops. 
The change in the central entropy values would only affect 
the limit on the entropy of the gas in infalling substructures 
that is allowed to penetrate the core, without affecting the 
general development of the profiles that we have described 
by the analysis of our results. 

Gas cooling, star formation and feedback processes from 
supernovae explosions associated to galaxy formation are 
taken into account in the semi-analytic model. The lack of 
AGN feedback in our model oblige us to use simple rules for 
the quenching of star formation to avoid overcooling in the 
semi-analytical prescription, and therefore we do not have a 
robust control of the cooled gas fraction. The impact of AGN 
feedback on galaxy properties is currently being investigated 
(Lagos, Cora & Padilla, in prep.); this improved semi-anlytic 
model will be used to study the influence of AGN feedback 
on the ICM metal enrichment. 

In the case of cool-core clusters, the formation and evo- 
lution of the central brigh t est galaxy plays an im portant 
role dBohringer et al.l 12004 |Pe Grandi et al.1 120041 '). Its ef- 
fect will be properly modelled with the inclusion of AGN 
feedback that mainly affects star formation in masive halos; 
this will allow to evaluate its influence on the development of 
the ICM metallicity, and to quantify the relative importance 
with respect to gas dynamical effects. 

The semi-analytic model does not consider tidal dis- 
ruption of cluster galaxies which give raise to intracluster 
stars, as shown by high resolution iV-body/SPH simulations 
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llMurante et al.l |2004| ; ISommer-Larsen. Romeo fe Portinaril 
l2005h . In our model, the potential intracluster stars re- 
main within the galaxies; their ejected supernovae prod- 
ucts are affected by the feedback process considered 
instead of being injected directly into the intracluster 
gas, the latter being the main advantange of the un- 
bound component in the enrichment process of the ICM 
l|Zaritskv, Gonzalez fe Zabludofj |2004| ). Thus, we can con- 
sider that our model underestimates the contribution of po- 
tential intracluster stars to the metal budget of the ICM. On 
the other hand, numerical simulations show that unbound 
stars (at least ~ 10 per cent of the cluster stellar mass) 
accumulate within the cluster generating intracluster light 
with a shallower radial profile than that of the bound com- 
ponent within ~ 0.3-Rvir! besides, they are on average older 
than the stars in cluster galaxies, which iron ejection rates 
already peak at high redshift. Hence, based on these results, 
we do not expect that intracluster stars drive a significant 
evolution of the ICM iron abundance profile at low redshifts. 

Ram pressure stripping is another process that might 
change the shape of the slope of the abundance pattern in 
the innermost part of the ICM. Its effect is present on the 
least bound gas particles in the hydrodynamical simulations. 
However, our semi-analytic model does not currently include 
the ram pressure stripping on the cold gas of cluster galaxies, 
which would affect the cold gas reservoir available for star 
formation and the timing of galactic mass loss. The evalua- 
tion of the impact of this mechanism will be part of future 
work. 
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